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ABSTRACT 


This paper presents a rapid procedure for generating c-type cascade grids 
suitable for viscous flow computations in turbomachinery blade rows. The 
resulting mesh is periodic from one blade passage to the next, nearly 
orthogonal, and continuous across the wake downstream of a blade. The 
procedure employs a pair of conformal mappings that take the exterior of the 
cascade into the interior of an infinite strip with curved boundaries. The 
final transformation to a rectangular computational domain is accomplished 
numerically. The boundary values are obtained from a panel solution of an 
integral equation and the interior values by a rapid ADI solution of Laplace’s 
equation. Examples of C-type grids are presented for both compressor and 
turbine blades and the extension of the procedure to three dimensions is 
briefly outlined. 
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Most of the coordinate systems in current use for turbomachinery flow 
computations are of one of three types. The channel gria has one family of 
lines starting upstream, passing through the blade rows, and continuing on 
downstream. The O-type grid has one family of lines that form closed loops 
around the blades. Finally, the C-type grid has one family of lines that wrap 
around the blade leading edge and continue on downstream. While the channel 
grid can be aligned with the flow and is fairly easy to generate, the 
resolution around the leading edge is usually poor and a choice usually has to 
be made between periodicity or near orthogonality for highly staggered 
cascades. Although the O-type grid provides excellent resolution around 
leading and trailing edges and may be both periodic and orthogonal, in general 
there is no mesh line aligned with the downstream flow and, hence, it is 
unsuitable for viscous computation. The C-type grid, on the other hand, 
appears to be a good choice for viscous flows. It provides good leading edge 
resolution, it can be both periodic arri orthogonal, and it can be aligned with 
the downstream flow. This paper presents a rapid procedure for generating 
such C-type grids. 
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The procedure starts with a conformal mapping which takes the exterior of a 
cascade of semi-infinite flat plates in the Z-plane into the interior of the 
unit circle in the W-plane. upstream infinity raps to the origin and 
downstream infinity to +1 on the real axis. This mapping is a limiting form 
of the standard rapping for a cascade of finite-chord flat plates (1) . When 
this mapping is applied to a real geometry, such as the turbine cascade in the 
figure, the semi-infinite flat plate is taken to run from a point 
inside the leading edge through the downstream end of the wake. 
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The mapping of the turbine cascade and wakes of the preceding figure produces 

the highly distorted "circle" in the adjacent figure. Note that the contour 

actually crosses the real axis twice between zero and one. The next napping 

takes the interior of the unit circle in the W plane, with a branch cut from 

zero to one along the real axis, to the interior of the infinite strip between 

the real axis and ( -i j ) in the ? -plane. The upper and lower sides of 

the wake at downstream infinity are mapped to plus and minus infinity, 

respectively, while upstream infinity maps to the origin. Since W is a 
2 

function of C , reflection of C through the origin leaves W unchanged. 

CONTOUR IN W PLANE 
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The image of the cascade of turbine blades and wakes in the C -plane is a 
pair of parallel straight lines coiinected by a roughly S-shaped curve. In 
actual practice W is eliminated between the two functions and the transfor- 
mation from Z to C is obtained by Newton iteration proceeding from point to 
point around the contour. To insure that the branch cuts of the logarithms 
are never crossed, the imaginary parts of these logarithms are saved. 

Whenever the change in either of these quantities between adjacent points 
exceeds ± it , the computed value of the logarithm at the new point is 
incremented by +2 mi, i.e. in the opposite direction. 


CONTOUR IN o PLANE 


The final mapping transforms the infinite strip in the ? -plane, bounded Vm 
by the blade-wake contour and its reflection through the origin, into a 
rectangular domain with coordinates F = £ + in* If we let F be the 
complex potential for flow through the strip and require F( c ) = -F(- C ) 
together with n = -1 along the contour, then we can write F as a contour 
integral. Here C, 6 , and h are, respectively, the complex velocity, flow 
angle, and normal channel width in the far fiela. The figure shows the 
formation of an integral equation for the unknown vortex density q . The 
source density q^ is set to cancel the normal component of the velocity 
C. Here s is distance along contour. A simple panel method, with flat 
panels and locally constant q fc and q n , is used to solve for and 

then to find £ as a function of Z along the contour. 

FORMATION OF INTE6RAL EQUATION ^ 

COMPLEX POTENTIAL 
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Generation of the grid in the rectangular E , n space proceeds in two 
stages. First points are located on the boundaries such that the grid in the 
cascade plane is periodic and continuous across the wake. Periodicity is 
enforced by distributing points symmetrically about the origin along the 
£ -axis. Continuity across the wake, away from the trailing edge, is 
achieved by selecting a constant mesh size A E, for this region such that the 
spacing in the cascade plane is an integer fraction of the staggered distance 
s sin | 3 W | , where s is the blade pitch and B w is the wake angle. The 
values of Z along the boundary are then found by inverse interpolation in 
the solution of E, vs Z .In order to enforce continuity near the 
trailing edge, a local straining is first introduced that places a point 
directly at the trailing edge. Then pairs of neighboring points across the 
wake are adjusted until their images in the cascade plane coincide. The 
distribution of points with n at the two ends of the region is arbitrary. 
Uniform spacing can be used for inviscid flows, and boundary layer stretching 
can be used to cluster points near the blade surface and wake for viscous 
flows. Once the boundary values of Z are specified, interior values are 
found by solving the complex Laplace equation by a cyclic ADI relaxation 
scheme which has the symmetry properties of Z built in. Estimates of the 
maximum and minimum eigenvalues of the matrix are used to obtain near optimum 
values of the acceleration parameters (2) . Convergence to the round-off error 
limit with seven place arithmetic is typically obtained in six to twelve 
iterations, even for cases where the maximum and minimum eigenvalues differ by 
five orders of magnitude. 


GRID GENERATION 
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ARBITRARY DISTRIBUTION ALONG n 
I INTERIOR VALUES IN ( { , n ) SPACE 


SOLVED BY CYCLIC ADI RELAXATION 


This figure shows the grid distribution in the z, -plane. Note that the 
upper boundary in the plot, which is found by the symmetric ADI solution of 
Laplace’s equation, naps into the upper and lower periodic lines in the 
cascade plane. 


GRID IN c PLANE 




The final grid in the cascade plane is obtained by conformal mapping of the 


solution in the Z -plane using the two analytic functions previously 
introduced. This figure shows the grid distribution for the cascade of 
turbine blades. Note that the continuity across the wake was obtained at the 
expense of a small amount of nonorthogonality. The rounded cap at the 
upstream boundary was obtained by extrapolation from the next two inner loops. 
Generation of this grid (99 x 7 points) required about 1.4 sec. of CPU time on 
an IBM 3033 computer. 

GRID IN CASCADE PLANE 
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The last figure presents C grids for a compressor stator and a turbine rotor. 
The stator was designed to turn the flow to the axial direction, hence, there 
is zero stagger in the downstream boundary. The turbine rotor is a 
particularly difficult case as it was designed to produce 126 degrees of 
turning. In this case the imposition of continuity across the wake resulted 
in a significant change in slope. 

The extension of this procedure to the generation of three-dimensional 
turbomachinery grids should be relatively straightforward. First the spanwise 
direction is discretized by a number of coaxial, axisynmetric surfaces. Next, 
ana most difficult, the intersection of the blade with each of these surfaces 
is obtained in meridional ( m ) and tangential ( 0 ) coordinates. Since the 
geometry is periodic in 6 , these ( m , 6 ) coordinates can be fed into the 
current program to generate a C-grid on each of the axi-symmetric surfaces. 

For O-grids this has already been done by Dulikravich (3) . 
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